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The non-perturbative Csw determined by the Schrodinger functional (SF) method with the RG-improved gauge 
action in dynamical Nf — 3 QCD shows a finite volume eff'ect when the numerical simulations are carried out at 
a constant lattice size L/a. We remove the unwanted finite volume effect by keeping physical lattice extent L at 
a constant. The details of the method and the result obtained for non-perturbative Csw with a constant L are 
reported. 
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1. Introduction 

Last year we reported the existence of a sizable 
finite volume effect in the non-perturbative Csw 
determined by the Schrodinger functional (SF) 
method at a fixed lattice size (L/a=const) when 
the RG-improved gauge action is adopted in full 
three-flavor QCD This effect causes a con- 
stant deviation from the true non-perturbative 
value of Cgw, and yields 0(a) errors in physical 
observables. 

This year we explain our strategy for removing 
the unwanted finite volume effect from Csw and 
report the bare coupling dependence of the non- 
perturbative Csw for three-flavor dynamical QCD. 

2. Method and Simulations 

The basic strategy to obtain the non- 
perturbative Csw with the SF setup is described 
in ig. We used Eq. (2.17) of 01 for the SF bound- 
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ary coupling of the RG-improved gauge action. 
The polynomial Hybrid Monte Carlo (PHMC) al- 
gorithm j3j is employed for the numerical simula- 
tion of the three-flavor dynamical QCD. 

To remove the finite lattice size effect we fix the 
physical lattice extent L so that a/L corrections 
vanish with a. This can be achieved by the fol- 
lowing steps. 

(1) Define the reference lattice size L* /a at a bare 
coupling value gg . This also defines the reference 
physical lattice extent L* . 

(2) For other 5o's? t^i^ lattice sizes L* /a at each 
(7q are estimated by the beta function. 

(3) Determine Csw and Kc via the PCAC condition 
with the SF setup at each gQ on the corresponding 
lattice size L* /a. 

In actual numerical simulations we employ 
the universal two-loop beta function to estimate 
L* /a. This causes a systematic error in the non- 
perturbative Csw, which could be large in the 
strong coupling region. This systematic error, 
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however, can be partly avoided by defining L* at 
the largest coupling value. We can also examine 
the magnitude of the error by investigating the 
L/a dependences because the error only appears 
through the incorrect value of L* /a. 

With the above consideration we take the fol- 
lowing setup in the simulations. 

(i) We define L* so that L*/a — 6 holds at 
/3 = 6/g2 = 1.9. 

(ii) We employ both Lja = 8 and 6 at = 2.0, 
and L/a = % for other weaker couplings. 

(iii) We correct Csw and at L/a to those &l L* /a 
using the one-loop correction calculated with the 
SF setup in finite lattice volume 0]. 
Although the perturbative error from setup (iii) 
appearers in strong coupling region, the magni- 
tude can be reduced when L* /a ~ L/a holds and 
L/a dependence is small. The small L/a depen- 
dence is confirmed at /? = 2.0 as shown later. We 
will briefly discuss the possible resolution of the 
systematic error from the use of perturbative beta 
function in the last section. 

We use the following non-perturbative tuning 
condition for and k. 



M{gl,L/a,n, 
AM{g^, L/a, 



= 0, 

) = o. 



(1) 



where M{gQ, L/a, k, Csw) is the bare PCAC quark 
mass and AM(gg, L/a, k, Cgw) is the difference of 
the two bare PCAC quark masses with different 
definition. The actual definitions are explained 
in [2]. The difference from the tuning condition 
we employed last year is the absence of the tree- 
level finite PCAC quark mass which arises from 
the SF boundary for a finite lattice extent. The 
tree-level finite volume effect is tuned to the cor- 
responding value at L*/a by setup (iii) together 
with the one-loop finite volume effect. 

As shown in setup (iii), we use the following 
one-loop perturbative correction to convert Cgw 
and K at L/a to those at L*/a at each gg; 



Kc{L* / a) = Kc{L/a) + SudL/a; L* / a), 



(2) 



Csw{L*/a) = Csw{L/a) + Scsw{L/a; L* /a) 
5k,c{L/ a;L* / a) and Scsw{L/a;L* /a) are defined as 
SKc{L/a;L*/a)^[Kc{L*/a)-Kc{L/a)]^^_^^^ , 



f„,= 1.50 
cw=1.60 
c™=1.70 
c™=1.80 
c™=l|.90 



(5csw(i/a;i*/a) = [csw(i*/a)-Csw(i/a)]i_, 



(3) 
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Figure 1. M vs AM at /3 = 2.0, L/a = 8. 

where the right hand side of Eq. ^ is calculated 
up to one-loop level. 

3. Results 

We have carried out the numerical simulation 
in the region /? = 12.0 — 1.9 with the setup ex- 
plained in Sec. 121 It was found that below (3 = 2.2 
the PHMC algorithm failed at vanishing or neg- 
ative PCAC quark masses due to large quantum 
fluctuations. We extrapolated them from M > Q 
to M = as shown in Fig. It was also found 
that below (3 = 1.9 with L/a = 8 it became dif- 
ficult to satisfy the condition AM — 0, which 
is why we employ L*/a = 6 at /3 = 1.9 for the 
definition of the reference physical lattice extent. 

Figures 121 and 13 are the result for Cgw and k^- 
Open circles are the result without L/a matching. 
Triangles are after tuning to L* /a by the one-loop 
correction (or interpolation at gQ — 3.0). For Csw 
we observe large one-loop corrections in the re- 
gion 1.0 < (7o < 2.7. This is due to rather strong 
coupling and a large discrepancy between L/a 
and L* /a. We discard the data from the func- 
tional fit. At strong couplings, g^ — 6/2.2 and 
6/2.1, the one-loop correction is rather small since 
L/a ^ L*/a. The L/a dependence at gQ — 3.0 
(/3 = 2.0) is negligible compared to the statis- 
tical error (upper circle: L/a = 8, lower circle: 
L/a = 6). We have 6 < L*/a < 8 via the two- 
loop beta-function at g^ = 3.0, we simply inter- 
polate Csw to L*/a using both lattice size data. 
While uncorrected data undershoots the one-loop 
result for L/a = oo (dotted line) in the weak cou- 
pling region those with one-loop correction 
gradually approaches it {L* /a is essentially infi- 
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uncorrected. 

1-loop corr. ((3=1.9,L*/o=6) (2-loop ^-func.) ' ^ 

1 -loop corr. 03=1.9.L*/o=6) (2-loop ^-func.) ' ^ 

1-loop <:„ at Mj=-> 
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Figure 2. gig dependence of Csw Open circles 
are the result before L/a tuning. Open and filled 
triangles are tuned to L* /a. Filled triangles are 
used for curve fitting. Solid line is the fit result. 
Dotted line shows the one-loop result at L/a=oo. 
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Figure 3. Same as Fig. |21 but for Kc- 

nite in the region). For Kc the one-loop correc- 
tion is rather small. We conclude that the non- 
perturbative results in the weak coupling region 
is consistent with the one-loop one for both Csw 
and Kc- We fit the filled triangles to obtain the 
functional form in for Csw and k^- 

After the curve fitting including the knowledge 
of the one-loop perturbative coefficient we obtain 

Csw(5o) = 1 + 0-11350 + 0.0209(72)54 
-fG.0047(27)ff6, 



(4) 



with xVDOF = 
Keigl) = 0.125 



0.5754, and 
f (3.681192 xlO 



(5) 



+ (0.141(43) X 10 
+ (0.124(67) X 10 
- (0.049(21) X 10 

with x^/DOF = 1.0014. All data are statistically 
independent and x^/DOF results in reasonable 
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values. We show statistical errors only. The sys- 
tematic error from the use of perturbative beta 
function is not estimated. We, however, consider 
that the systematic error is minimized since the 
slope in = 3.0 — 3.15 of Csw should be cor- 
rect because of the negligible L/a dependence at 
gl = 3.0. 

4. Summary 

We have extracted the non-perturbative Csw as 
a function of (?q for Nj — 3 with the RG-improved 
gauge action. The functional form is shown in 
Eq. Q). In deriving the non-perturbative values 
at a fixed physical lattice extent L*, we used the 
two-loop beta-function to obtain the correspond- 
ing lattice size L*/a at a given g^. This may 
contradict the word "non-perturbative" . In or- 
der to improve this point we should iteratively 
execute our procedure with the non-perturbative 
beta-function that is measured with the non- 
perturbative Csw extracted from the previous it- 
eration. We leave this iterative study for future 
studies. The first step to obtain Nf = 2+1 
hadron mass spectrum has been started using the 
non-perturbative Csw of Eq. I^J ^ . 
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